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ABSTRACT 

CZ i 

We present two-dimensional (axisymmetric) neutrino-hydrodynamic simulations of the long-time accretion phase of a 15 M Q progen- 
itor star after core bounce and before the launch of a supernova explosion, when non-radial hydrodynamic instabilities like convection 
occur in different regions of the collapsing stellar core and the standing accretion shock instability (SASI) leads to large-amplitude 
oscillations of the stalled shock with a period of tens of milliseconds. Our simulations were performed with the Prometheus- Vertex 
code, which includes a multi-flavor, energy-dependent neutrino transport scheme and employs an effective relativistic gravitational 
potential. Testing the influence of a stiff and a soft equation of state for hot neutron star matter, we find that the non-radial mass 
motions in the supernova core impose a time variability on the neutrino and gravitational-wave signals with larger amplitudes, as 
well as higher frequencies in the case of a more compact nascent neutron star. After the prompt shock-breakout burst of electron 
neutrinos, a more compact accreting remnant produces higher neutrino luminosities and higher mean neutrino energies. The observ- 
' able neutrino emission in the SASI sloshing direction exhibits a modulation of several ten percent in the luminosities and around 

1 MeV in the mean energies with most power at typical SASI frequencies between roughly 20 and 100 Hz. The modulation is caused 
by quasi-periodic variations in the mass accretion rate of the neutron star in each hemisphere. At times later than ~50-100ms after 
bounce, the gravitational-wave amplitude is dominated by the growing low-frequency (S200 Hz) signal associated with anisotropic 
neutrino emission. A high-frequency wave signal results from nonradial gas flows in the outer layers of the anisotropically accreting 
neutron star. Right after bounce such nonradial mass motions occur due to prompt post-shock convection in both considered cases 
and contribute mostly to the early wave production around 100 Hz. Later they are instigated by the SASI and by convective overturn 
CO . that vigorously stir the neutrino-heating and cooling layers, and also by convective activity developing below the neutrinosphere. The 

gravitational-wave power then peaks at about 300-800 Hz, connected to changes in the mass quadrupole moment on a timescale of 
milliseconds. Distinctively higher spectral frequencies originate from the more compact and more rapidly contracting neutron star. 
Both the neutrino and gravitational-wave emission therefore carry information that is characteristic of the properties of the nuclear 
equation of state in the hot remnant. The detectability of the SASI effects in the neutrino and gravitational-wave signals is briefly 
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, 1 . Introduction that deform and stir the condensing central compact remnant of 

jrt ■ the explosion. 

- - i Neutrinos and gravitational waves are the most direct potential 

probes of the processes that occur deep inside of a dying star, The P r °mpt burst of electron neutrinos, for example, signals 
accompanying or causing the initiation of the stellar explosion. the breakout of the supernova shock from the neutrinosphere and 
Neutrinos were already detected in connection with supernova the shock heating of matter in this region. The neutrino emission 
SN 1987A (Bionta et al. 1987, Hirata et al. 1987, Alexeyev et after core bounce is a sensitive probe of the mass accretion rate 
al. 1988), although with poor statistics so that the extraction of on the forming neutron star, hence of the density structure of the 
information for constraining the explosion mechanism was not falling layers of the dying star (e.g., Liebendorfer et al. 2003, 
possible. The simultaneous measurement of signals of both types Buras et al. 2006b). The onset of the explosion is expected to 
remains a very realistic hope for the next Galactic supernova. show U P as a more °r less sudden drop in the mass accretion rate 
Capturing neutrinos and gravitational waves from the same of the formin g neutron star - A soft supernuclear equation of state 
source has the advantage of providing complementary insight wl11 lead to a more com P act and hotter remnant, radiating higher 
into the conditions of the stellar core. While neutrino signals re- neutrmo luminosities and more energetic neutrinos (e.g., Janka 
fleet the density structure and thermodynamic conditions in the et aL 2005 ' Marek 2007 )- And a P osslble P hase transition to non- 
high-density plasma of the collapsing core and forming neutron nucleomc matter in the neutron star core may impose character- 
star, gravitational waves carry crucial information about the dy- lstlc features llke a second neutnno burst ( Sa 8 ert et aL 2008 ) or ' 
namics and nonradial motions of the stellar matter, e.g. of its ro- lf triggering the collapse to a black hole, may cause an abrupt 
tational state or of hydrodynamic instabilities such as convection termination of the neutrino emission (e.g., Burrows 1988; Keil 
3 3 & Janka 1995; Sumiyoshi et al. 2006, 2007; Fischer et al. 2008). 

Send offprint requests to: H.-Th. Janka, e-mail: Theoretical work on the gravitational-wave signals from stel- 

thj@mpa-garching.mpg.de lar core collapse and explosion has a long history of successively 
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refined numerical models. In particular the infall and bounce 
phases, which are theoretically relatively well understood parts 
of the evolution, have received a lot of interest, because a strong 
and characteristic signature could make them a promising source 
of a detectable gravitational-wave burst (for a review-like intro- 
duction to the subject and a nearly complete list of publications, 
see Dimmelmeier et al. 2008). For this to be the case, the core 
of the progenitor star must develop a sufficiently large deforma- 
tion during its infall, and for that it must rotate enough rapidly. 
This, however, does not seem to be compatible with predictions 
from the latest generation of stellar evolution models for the vast 
majority of massive stars, which are expected to have lost most 
of their angular momentum before collapse (Heger, Woosley, & 
Spruit 2005). Only in rare, very special cases, possibly account- 
ing for the few tenths of a percent of all stellar core collapses 
that produce gamma-ray bursts, such stars seem to be able to 
retain a high angular momentum in their core and to thus pro- 
duce relativistic jets and highly asymmetric and extremely en- 
ergetic, probably magnetohydrodynamically driven explosions 
(for a recent review, see e.g. Woosley & Bloom 2006). If this 
hypothetical connection was true, the core bounce phase would 
not really offer grand perspectives for the measurement of grav- 
itational waves. 

In contrast, only relatively little work has so far been done 
on determining the wave signals from the post-bounce evolu- 
tion of a supernova, although these signals are likely to carry 
important information about the still incompletely understood 
explosion mechanism and the associated core dynamics (see the 
review by Ott 2008). Miiller & Janka (1997) and later Miiller 
et al. (2004) on the basis of significantly improved numerical 
models, showed that in the case of delayed, neutrino-driven ex- 
plosions convective overturn behind the stalled shock as well as 
convection inside the nascent neutron star can account for siz- 
able gravitational-wave signals, which should be detectable with 
a high probability from a Galactic supernova when the Advanced 
Laser Interferometer Gravitational-Wave Observatory (LIGO II) 
is running (see also Fryer et al. 2002, 2004). 

While the delayed neutrino-heating mechanism relies on the 
support by strong nonradial hydrodynamic instabilities in the re- 
gion of neutrino-energy deposition, magnetohydrodynamic ex- 
plosions, if linked to rapid rotation, are expected to develop rela- 
tively soon after core bounce and thus to occur faster than the 
growth of the mentioned nonradial instabilities. They tap the 
reservoir of differential rotation in a collapsing stellar environ- 
ment and therefore require rapidly spinning cores (see Burrows 
et al. 2007b, Thompson et al. 2005). Thus they are good can- 
didates for sizable gravitational-wave pulses from the moment 
of core bounce (e.g., Kotake et al. 2006, Ott et al. 2004). The 
newly proposed acoustic explosion mechanism, in which large- 
amplitude core gravity-mode oscillations of the neutron star con- 
vert accretion power to pressure and shock waves that feed the 
supernova shock with acoustic energy (Burrows et al. 2006, 
2007a), is predicted to be associated with enormous and very 
characteristic gravitational-wave activity due to the fast periodic 
movement of roughly a solar mass of dense matter at late times 
(~1 second) after core bounce (Ott et al. 2006). 

In the present paper our focus is on an analysis of the fea- 
tures in the neutrino and gravitational-wave signals that might 
provide evidence for the action of the so-called standing accre- 
tion shock instability (SASI), which has been shown to lead to 
low-^ mode (in terms of an expansion in spherical harmonics 
with order {), large-amplitude nonradial shock deformation and 
violent sloshing motions of the stalled supernova shock (Blondin 
et al. 2003; Blondin & Mezzacappa 2007; Bruenn et al. 2006; 



Scheck et al. 2004, 2006, 2008; Ohnishi et al. 2006, Foglizzo 
et al. 2007; Yamasaki & Foglizzo 2008). This SASI activity is 
found to reach the nonlinear regime at roughly 100 ms after core 
bounce and to grow in amplitude over possibly hundreds of mil- 
liseconds (Scheck et al. 2008; Marek & Janka 2007; Burrows 
et al. 2006, 2007a). It does not only act as seed of powerful 
secondary convection but can also provide crucial aid for the 
neutrino-heating mechanism by pushing the accretion shock to 
larger radii and by thus stretching the time accreted matter is ex- 
posed to neutrino heating in the gain layer (Buras et al. 2006b, 
Scheck et al. 2008, Marek & Janka 2007, Murphy & Burrows 
2008). Moreover, the SASI is found to play the driving force 
of the g-modes pulsations of the neutron star core that are the 
essential ingredient of the acoustic mechanism. 

Observational signatures of the presence of the SASI would 
therefore be extremely important for our understanding of how 
massive stars begin their explosion. Supernova asymmetries and 
pulsar kicks are one, yet not unambiguous observational hint. 
Neutrinos and gravitational waves, which originate directly from 
the region where the blast is initiated, however, may remain the 
only way to obtain direct information. It is therefore a highly 
relevant question to ask whether any characteristic structures are 
imprinted on the neutrino and gravitational-wave emission by 
the SASI activity in the supernova core. In contrast to convec- 
tion, which exhibits the fastest growth for the higher-^ modes 
(see Foglizzo et al. 2006), the SASI is expected to possess the 
largest growth rates for the dipolar and quadrupolar deforma- 
tions (corresponding to £ = 1,2; Blondin & Mezzacappa 2006; 
Foglizzo et al. 2007; Yamasaki & Foglizzo 2008; Foglizzo 2001, 
2002). Even in the fully nonlinear situation the geometry and 
motion of the shock and post-shock layer are found to be gov- 
erned by these lowest modes. 

We present here an analysis of the neutrino and gravitational- 
wave signals that are calculated on the basis of the two- 
dimensional post-bounce and pre-explosion simulations of a 
15 M star recently published by Marek & Janka (2007). We 
constrain ourselves to two nonrotating models, which allow us to 
discuss the differences that can be expected from a stiff and a soft 
nuclear equation of state (results of the corresponding ID simu- 
lations can also be found in Marek & Janka 2007). We find that 
the SASI sloshing of the shock and the associated quasi-periodic 
mass motions lead to a time-modulation of the neutrino emission 
and to gravitational-wave amplitudes whose size and character- 
istic frequency depend on the compactness of the proto-neutron 
star during the first half of a second after core bounce. Though 
the SASI contributions to the power spectra are superimposed 
by a significant high-frequency "noise" due to convective fluctu- 
ations, at least the combination of measurements should make it 
possible to identify the SASI activity in the supernova core. 

Our paper is organized as follows. In Sect. [2] we will 
briefly outline the main numerical and physics ingredients of 
our simulations, and the basic features of the two simulations 
we compare. In Sect. [3] we discuss our results with respect 
to the SASI effects on the shock motion, neutrino emission, 
and gravitational-wave signal, and in Sect. [4] we will summa- 
rize our findings and draw conclusions, including a discussion 
whether the SASI modulations are detectable by neutrino and 
gravitational-wave experiments. 

2. Code, input, models 

The 2D simulations discussed in this paper were performed with 
the Prometheus-Vertex code, whose numerical aspects and the 
implemented microphysics were described by Rampp & Janka 
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(2002) and Buras et al. (2006a), and the publications quoted 
in those papers. The detailed list of ingredients was also pro- 
vided by Marek & Janka (2007), where the simulations were 
already introduced as Models M15LS-2D and M15HW-2D and 
compared with other cases in a greater set of calculations. We 
therefore repeat only a few essential aspects of immediate rele- 
vance here and refer the reader to Sect. 2 of the Marek & Janka 
(2007) paper for more complete information. 

The hydrodynamics module of the code is based on a con- 
servative and explicit Eulerian implementation of a Godunov- 
type scheme with higher order spatial and temporal accuracy. 
It solves the nonrelativistic equations of motions for the stellar 
fluid, whose self-gravity is described by an "effective relativistic 
potential" for an approximative treatment of general relativis- 
tic gravity (see Marek et al. 2006; the discussed simulations 
were performed with the potential of Case A from this work). 
The neutrino transport, which is coupled to the hydrodynamics 
part via lepton number, energy, and momentum source terms, is 
computed with our "ray-by-ray plus scheme" (see Buras et al. 
2006a). It accounts for the full neutrino-energy dependence of 
the transport but treats its dependence on the direction of the neu- 
trino momentum in an approximative way, which is numerically 
less demanding and more efficient than a full multi-dimensional 
version of the transport. Recent multi-angle simulations in two- 
dimensional situations (Ott et al. 2008) — though done at the 
expense of a sophisticated description of the energy dependence 
— show that a detailed angular treatment produces considerably 
less lateral smearing of the outward directed radiation field than 
flux-limited diffusion. We suspect that our ray-by-ray descrip- 
tion compares much better with multi-angle results than flux- 
limited diffusion does. 

The progenitor star used for our simulations was model 
sl5s7b2 from Woosley & Weaver (1995), which is a standard 
nonrotating 15 M Q star widely used for supernova simulations. 
We employed two different nuclear equations of state (EoS) for 
our studies: (1) a soft version of the Lattimer & Swesty (1991) 
EoS ("L&S EoS") with an incompressibility modulus of bulk 
nuclear matter of 180MeV and a symmetry energy parameter 
of 29.3 MeV, and (2) the considerably stiffer EoS of Hillebrandt 
& Wolff (1985; "H&W EoS"; see also Hillebrandt et al. 1984), 
whose parameter values are 263 MeV and 32.9 MeV, respec- 
tivel)o3- The former leads to a radius of about 12 km for cold 
neutron stars with a "typical" (gravitational) mass of 1.4 M , 
whereas this radius is roughly 14 km in the second case. The 
Hillebrandt & Wolff EoS is based on a largely different model- 
ing approach for inhomogeneous nucleon matter than the L&S 
EoS and the more recent EoS of Shen et al. (1998). It employs 
a Hartree-Fock calculation in contrast to the compressible liq- 
uid drop model of the L&S EoS and the relativistic mean field 



1 Note that below a certain density, which is typically chosen to be 
10" gcnT 3 after core bounce, we replace the high-density EoSs by an 
ideal-gas equation of state with electrons, positrons, photons, and a mix- 
ture of classical, nonrelativistic Boltzmann gases for nucleons, alpha 
particles, and 14 kinds of heavier nuclei (see Marek & Janka 2007 for 
more information). 

2 In order to save computer time, in particular during the simulation 
phase right after core bounce when the timesteps are constrained to very 
low values, the model with the H&W EoS was performed with the as- 
sumption of equatorial symmetry until 125.3 ms post bounce. Only af- 
terwards it was continued with a full 180° grid. Tests showed that this 
had no important influence on the results, neither for the evolution of the 
shock radius, nor for the growth and development of hydrodynamic in- 
stabilities. In particular, we did not find any qualitative difference in the 
SASI modes, also for odd values of I, at later times (see Marek 2007). 



description of the Shen et al. EoS. The three EoSs yield sig- 
nificantly different results in ID core-collapse simulations with 
respect to the shock formation point, the luminosities and mean 
energies of the radiated neutrinos, and the evolution of the shock 
radius and neutron star radius after bounce (see Janka et al. 
2005; Marek 2007; Figs. 6 and 7 in Janka et al. 2007). Since 
the EoS of Shen et al. yields intermediate values for many of 
these quantities, we consider the soft Lattimer & Swesty EoS on 
the one hand and the stiff Hillebrandt & Wolff EoS on the other 
as two cases that roughly span the range of extreme possibilities 
for baryonic matter around and above nuclear saturation density 
near core bounce and shortly afterwards. 

3. Results 

In the following we will compare our two simulations for these 
equations of state in detail, focussing on the observable signals 
whose properties depend on the asymmetries due to hydrody- 
namic instabilities in the supernova core after bounce. 

3.1. Hydrodynamic instabilities and shock motion 

Basically one can discriminate three regions and episodes 
of nonradial hydrodynamic instability in the stellar core af- 
ter bounce: (i) prompt post-shock convection (Epstein 1979, 
Burrows & Fryxell 1993, Janka & Miiller 1996, Miiller & Janka 
1997, Swesty & Myra 2005), (ii) convection inside the nascent 
neutron star, i.e. below the neutrinosphere of the electron neutri- 
nos (Burrows 1987, Keil et al. 1996, Buras et al. 2006b, Dessart 
et al. 2006), and (iii) convective overturn in the neutrino-heating 
layer between the gain radius and the stalled supernova shock 
(Herant et al. 1994; Burrows et al. 1995; Janka & Miiller 1996; 
Fryer & Warren 2002, 2004) as well as SASI activity (Blondin 
et al. 2003, Blondin & Mezzacappa 2006, Ohnishi et al. 2006, 
Scheck et al. 2008). These regions can be identified for both 2D 
simulations in Fig.Q] 

Region (i) arises when the newly formed supernova shock 
propagates outward through the infalling stellar iron core 
and experiences massive energy losses by nuclear photo- 
disintegrations. This weakens the shock so that it leaves be- 
hind a layer with negative entropy gradient (before it encoun- 
ters infalling material with even lower densities and therefore 
produces rising gas entropies again in spite of its further decel- 
eration). In addition, a Y e minimum is established near an en- 
closed mass of roughly 1 M when a large number of newly pro- 
duced electron neutrinos begins to stream away from the neutri- 
nosphere in a luminous shock-breakout burst. The electron frac- 
tion forms a trough around the radial position where the shock 
makes the transition from the neutrino-opaque to the neutrino- 
transparent regime (Fig. |2|. The electron fraction decreases to- 
wards this minimum, because electron-capture neutrinos escape 
more slowly from regions with higher densities, and because 
the conversion of electrons to neutrinos proceeds more slowly 
at lower densities. In both simulations the negative entropy gra- 
dient overlaps partly with the negative Y e gradient. 

The layer where both the entropy and the Y e gradients are 
negative is found to be Ledoux unstable (marked by the grey 
vertical bars in Fig. [2]). This region lies between the first local 
maximum and the following local minimum of the entropy pro- 
file. It contains considerably more mass in the model computed 
with the H&W EoS, where it is bounded by the mass coordi- 
nates of 0.63 M Q and 0.79 M . At higher enclosed masses, the 
entropy rises slowly out of its local minimum. In contrast, in the 
L&S case the boundaries are at 0.60 M and 0.71 M , and the 
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Fig. 1. Regions of mass motions due to hydrodynamic instabilities in the supernova core. The plots show color coded the lateral 
velocity (top) and the radial velocity (bottom) in the equatorial plane of the polar grid as functions of post-bounce time and radius 
for the 2D simulation with the L&S EoS (left) and for the simulation with the H&W EoS (right). The displayed range of velocity 
values is limited to +2 x 10 8 cms _1 . The white dashed lines mark (with decreasing radius) the locations of densities 10 10 , 10 12 , 
and 10 14 gcm -3 . Prompt post-shock convection is strongest at r ^ 15 km before t ~ 30 ms after core bounce, proto-neutron star 
convection occurs later at r ^ 30 km and p > 10 12 g crrT 3 , and SASI and convective activity behind the shock are visible at r k 30- 
40 km. 
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Fig. 2. Entropy (black) and electron fraction (red) profiles (the plotted values correspond to 10 x Y e ) as functions of enclosed mass 
at 10 ms after core bounce in spherically symmetric counterparts of the two 2D simulations discussed in this paper. The shock is the 
entropy jump near 1.1 M©. The vertical grey bars indicate the mass regions that are convectively unstable according to the Ledoux 
criterion. 



increase of the entropy in the overlying shells is much steeper. 
In the H&W model the thicker layer in mass corresponds also 
to a much wider radial domain at a larger distance from the stel- 
lar center. Because of higher growth rates the convective activity 
sets in more immediately after the shock passage and lateron be- 
comes stronger and encompasses a greater fraction of the stellar 
core. This can be seen at t S= 30 ms after bounce in Fig.Q] where 
the righthand panels for the H&W simulation in comparison to 
the lefthand panels show significantly more violent convective 



overturn. While the L&S model develops visible mass motions 
between ~15 km and ~40km, the H&W model shows stronger 
activity at radii of 15 km 5= r % 80 km. 

After some 10 ms the profiles of entropy and electron frac- 
tion are flattened and with the disappearance of the driving force 
the convective activity begins to calm down (this can be recog- 
nized better in the two lower panels of Fig. [TJ. Now, however, 
region (ii) is being formed. Neutrino transport begins to reduce 
the electron fraction deeper inside the nascent neutron star. As a 
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Fig. 3. Four representative snapshots from the 2D simulation with the L&S EoS at post-bounce times of 247 ms (top left), 255 ms ( top 
right), 322 ms (bottom left), and 375 ms (bottom right). The lefthand panel of each figure shows color-coded the entropy distribution, 
the righthand panel the radial velocity component with white and whitish hues denoting matter at or near rest; black arrows in the 
righthand panel indicate the direction of the velocity field in the post-shock region (arrows were plotted only in regions where the 
absolute values of the velocities were less than 2 x 10 9 cms _1 ). The vertical axis is the symmetry axis of the 2D simulation. The 
plots visualize the accretion funnels and expansion flows in the SASI layer, but the chosen color maps are unable to resolve the 
convective shell inside the nascent neutron star. 



consequence of this the lower boundary of the Y e trough and thus 
the negative Y e gradient moves gradually inward. At the same 
time the energy losses by neutrinos reduce the level of the en- 
tropy plateau at s ~ 5 fce/nucleon such that the plateau becomes 
wider and its inner edge also moves inward. These changes in 
the entropy and lepton number profiles in combination lead to 
the establishment (or/and maintenance) of a convectively un- 
stable layer inside the proto-neutron star. About 50-60 ms af- 
ter bounce, following a short, intermediate period of relative 
quiescence, both of our 2D models indeed show the reappear- 
ance of vivid overturn activity at densities above 10 12 gcm~ 3 
(Fig. [D- The convective transport of lepton number and entropy 
supports and enhances the mentioned structural changes so that 
a long-lasting zone with convective mixing develops (compare 
our Fig.Q]with Figs. 9 and 10 of Dessart et al. 2006). This zone 
encompasses a growing mass region within the broadening Y e 
trough. 

Region (iii) with nonradial hydrodynamic mass motions de- 
velops outside of the steep density gradient at the surface of 



the nascent neutron star, where the neutrinospheres of all neutri- 
nos are located. First signs of low-£ mode, low-amplitude SASI 
sloshing can be seen here already some 10 ms after bounce. At 
about 100 ms after bounce this activity strengthens significantly 
because at this time neutrino energy deposition in the gain layer 
has begun to build up a negative entropy gradient between the 
gain radius and the stagnant supernova shock. The neutrino- 
heating region is unstable to convection, and the presence of 
convective motions produces vorticity and entropy perturbations 
that feed back into the advective-acoustic cycle that is consid- 
ered as an explanation of the SASI phenomenon (see Foglizzo 
2001, 2002; Foglizzo et al. 2007). Thus the onset of convec- 
tion may lead to an amplification of the SASI growth. On the 
other hand, the violent shock oscillations that are characteristic 
of the fully developed SASI produce large entropy variations in 
the downstream region. These act as seeds for secondary convec- 
tion (Scheck et al. 2008) and have been found to aid the neutrino- 
heating mechanism for powering supernova explosions (Scheck 
et al. 2008, Marek & Janka 2007; see also Murphy & Burrows 
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Fig. 4. Left: Maximum shock radii (solid lines) and proto-neutron star radii (dashed lines) as functions of post-bounce time for the 
2D simulations with different nuclear equations of state. The neutron star radii are determined as the locations where the rest-mass 
density is equal to 10 1 1 g crrT 3 . Right: Shock contours at the different post-bounce times listed in the figure. The vertical axis of the 
plot is the symmetry axis of the simulation. 



2008). In the nonlinear regime convective and SASI activity are 
inseparably coupled, and a high-£ convective mode pattern oc- 
curs superimposed on the low-^ mode SASI deformations of the 
shock contour and of the post-shock region (for snapshots, see 
Fig. [3]and Marek & Janka 2007, Scheck et al. 2008). 

Although both simulations show the same basic features and 
evolutionary stages, the differences in details are interesting. 
As mentioned above, the H&W model develops more vigor- 
ous prompt post-shock convection in a region with higher mass 
and wider radial extension. In contrast, the SASI activity in this 
model is appreciably less strong than that in the calculation with 
the L&S EoS until roughly 350 ms after bounce. This can be seen 
in Fig. |4] where the maximum shock radius of the L&S case ex- 
hibits significantly bigger SASI amplitudes and correspondingly 
the shock contours show more extreme nonspherical deforma- 
tion (righthand panel of Fig. 0}. These differences grow during 
the nonlinear phase of the SASI between ~150 and 350 ms until 
at t <; 350 ms p.b. the SASI sloshing in the H&W model also 
gains more power and the maximum shock radii of both simula- 
tions become more similar again. 

We suspect that the more compact proto-neutron star for the 
softer L&S EoS leads to conditions that favor strong SASI activ- 
ity, possibly because of the more efficient neutrino heating (and 
more vigorous post-shock convection) as a consequence of the 
higher neutrino luminosities and harder neutrino spectra that are 
radiated from a more compact and hotter nascent neutron star 
(see Figs. [6] and [7] and Sect. |3.2| >. Another possible explanation 
in the context of the advective-acoustic cycle scenario may be 
different amplification factors of perturbations in the two mod- 
els. Due to the lack of good theoretical insight into the behavior 
of the SASI in the fully nonlinear regime, we do not see a way 
how to facilitate deeper understanding by further analysis. For 
both suggested explanations, however, one might expect that in 
the H&W simulation at late post-bounce times, when the proto- 
neutron star radius has contracted (Fig. |4]i and the neutrino- 
heating timescale has decreased (see Fig. 6 in Marek & Janka 
2007), the conditions for violent SASI activity have the tendency 
to improve. This would be consistent with our observed growth 
of the SASI amplitudes near the end of this simulation. 



In Fig. [5] lefthand panels, we plot the time-evolution of 
the amplitudes A\$, A2,o, and A3 0, of the decomposition of the 
angle-dependent shock radius R s (6, t) into spherical harmonics: 

00 

R s (0,t) = J]Ac S) (t)P° e (cos6) , (1) 

e=o 

where m — because of the axial symmetry of our models and 
P"(cos#) are the Legendre polynomials^. Figure [5] confirms our 
description above: the amplitudes for the L&S run are typically 
2-3 times larger than those of the H&W simulation. Only near 
the end of the computed evolution (at t <; 350 ms p.b.), A 10 and 
A3 of both models reach very similar values. This is not so for 
A-2,0, which exhibits oscillatory behavior on top of a clear slope 
to positive values in both simulations, signaling a trend in grow- 
ing prolate deformation of the shock. Until the end of the com- 
putational runs, however, this deformation remains considerably 
more pronounced in the case of the model with the L&S EoS. 

The righthand panels of Fig. [5] provide the Fourier trans- 
forms of the time-dependent mode amplitudes. One can see a 
broad peak of the Fourier spectra between about 10 and roughly 
100 Hz, followed by a steep decline towards higher frequencies. 
The t = 1 amplitude exhibits the clearest, jagged peak struc- 
ture at a frequency of 20-60 Hz for the L&S run, and a smoother 
peak at 30^0 Hz with a secondary one around 65 Hz for the 
H&W model. A 2,0 shows most power at frequencies around SO- 
SO Hz and 70-90 Hz in both simulations, and A3 possesses a 
broad Fourier maximum below ~100-130Hz, with an indica- 
tion of a peak at about 50 Hz in the H&W case. There is no 
really clear correlation of the SASI oscillation period with the 
compactness of the nascent neutron star or with the average ra- 
dius of the stalled shock except for, maybe, the slightly higher 

3 Note that we have renormalized our amplitudes A [0 compared to 
the coefficients ac used by Kotake et al. (2007) in order to directly dis- 
play in Fig.[5]the radius variations associated with the different spherical 
harmonics components. The coefficient Ao,o is identical with the angu- 
lar average of the shock radius: A .o(?) = (R s (t)) = (47r) _1 fdQR s (6,t), 
whose values and evolution are similar to those of the maximum shock 
radius plotted in Fig. [4] 
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mode towards the end of our simulations signals an increasing and persistent prolate deformation of the shock. Note that during 
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performed only with a 90° wedge from the pole to the equator, and we can therefore plot only the amplitudes with even t for this 
equatorially symmetric phase. Right: The Fourier transforms of the time-dependent mode amplitudes. 



frequency of the main I = 1 peak in the L&S run. Because of the shock region the appearance of well localized and very promi- 
time-dependent structure of the accretor and of the whole post- nent peaks in the spectrum might not be expected. In addition, 
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shock wobbling associated with convective mass motions pro- 
duces a background of short-wavelength "noise" at frequences 
up to roughly 200 Hz (Pig. [3j. Accordingly, in particular the am- 
plitudes for higher i values show power nearly evenly distributed 
over a wide range of frequencies. The Fourier transforms for 
^i,o> ^2,o> an d ^3,0 plotted in Fig. [5] nevertheless confirm the 
stronger SASI activity in the simulation with the L&S EoS: only 
for a few, narrow frequency windows the Fourier amplitudes of 
the H&W model become larger than those of the L&S run. 

In the following two sections we will turn to an analysis 
of the consequences of the described nonradial hydrodynamic 
instabilities, in particular also of the SASI, for observable sig- 
nals from the supernova core. One question of interest will be, 
whether the differences caused by the use of stiff or soft neutron 
star equations of state manifest themselves in some distinctive 
features of the signals. 

3.2. Neutrino emission 

The luminosities for neutrinos and antineutrinos of all flavors 
and for both 2D simulations, compared to the results of the corre- 
sponding ID models, are given as functions of post-bounce time 
in Fig.|6](muon and tau neutrinos and their antip articles, which 
we sometimes denote with v x , are treated in the same way, and 
the plots show the luminosities of one individual type of these 
heavy-lepton neutrinos). The luminosities displayed there are ac- 
tually not defined as the total rate of energy loss of the supernova 
core through neutrinos radiated in all directions. Instead, they 
represent the isotropic equivalent luminosities inferred by an ob- 
server at rest relative to the stellar center, viewing the source 
from the direction of the polar axis of the grid and from a large 
distance above the north pole (i.e., gravitational redshifting is 
taken into account). The results for an observer above the south 
pole look qualitatively and quantitatively very similar and lead 
to the same conclusions. 

One of the goals of the present work is an analysis of lumi- 
nosity fluctuations that are associated with modulations of the 
mass accretion rate of the nascent star due to SASI oscillations 
and convective overturn in the post-shock layer, which become 
very vigorous and create nonstationary conditions at later post- 
bounce times (see Sect. 13. II and Figs. [4] and [5]). Strong shock re- 
traction leads to a transient increase of the gas flow towards the 
neutron star and to the compression and enhanced cooling of the 
matter near the neutron star surface. In contrast, shock expansion 
has the opposite effect because it causes a deceleration of the 
infall or even outward acceleration of material that is accreted 
through the shock. Thus shock expansion stretches the time this 
matter remains in the gain layer and does not cool by neutrino 
emission (Marek & Janka 2007, Scheck et al. 2008, Murphy & 
Burrows 2008). Corresponding quasi-periodic fluctuations of the 
neutrino luminosities, strongest for electron neutrinos v e and an- 
tineutrinos v e (which are more abundantly produced in the layer 
of freshly accreted, hot matter near the neutron star surface) can 
be seen in all panels of Fig. [6] 

In the lefthand panels of Fig. [6] we have constructed the 
isotropic equivalent luminosities from the neutrino radiation 
leaving the stellar core in one angular grid bin near the north 
pole of the mesh, i.e., 

L v ,i(e M , f) = 4nr 2 F v (r, 6 j+ , , f) (2) 

for F v (r, 9j + i, f) being the energy-integrated neutrino energy flux 
in the rest frame at the cell-centered (pole-near) latitudinal angle 
6j + i of the polar grid of the simulation, at time t and a radius 



r that is large enough such that the luminosity is conserved at 
greater radial distances (for the numerical evaluation we have 
chosen r = 400 km). In contrast, for the righthand panels of 
this figure we have integrated the neutrino emission over the 
whole northern hemisphere and then rescaled the result to the 
full sphere, i.e., 

L v , 2 (f) = 2r 2 f dQF y (r,6,f) (3) 

■^hemisphere 

= 4nr 2 ]T (cos Q h i - cos 0j)F v (r, 6j_ i , t) , (4) 

7=1 

where the second line exploits axial symmetry and the sum is 
performed over the angular zones of one hemisphere. The former 
approach means that one assumes that an observer at a position 
near the polar axis receives the radiation emitted only from one 
pole-near lateral zone (in our treatment this means that the re- 
ceived radiation flux has the properties calculated in exactly one 
angular grid bin or angular "ray"), whereas the second evaluation 
implies the assumption that each unit of area on the hemisphere 
oriented towards the observer contributes to the observable lu- 
minosity with a weight defined by the (radial) flux calculated 
for the corresponding angular bin ("ray") of the polar grid in 
our ray-by-ray transport treatment. This assumption would be 
correct, of course, if a spherical source were radiating uniformly 
and isotropically in all directions, but in general the surface parts 
oriented with smaller angles to the observer direction contribute 
more strongly. For a sphere with a sharp radiating surface (in 
contrast to a radially extended neutrino-decoupling layer) that 
emits neutrinos with a (locally) isotropic intensity towards an 
observer at great distance on the axis of symmetry of the source, 
Eq. (O would contain a factor (2 cos 0) in the integrancfl If the 
intensity was a function of the latitudinal angle 6 of the stellar 
grid, for example, this factor would enhance the relative contri- 
butions of surface parts that are closer to the observer. Because 
of limb darkening the importance of the observer-near regions of 
the radiating neutrinosphere is further enhanced. 

Not solving the full multi-angle problem of neutrino trans- 
port in the 2D geometry but using our ray-by-ray transport ap- 
proximation does not allow us to exactly compute the direction 
dependence of the neutrino intensity radiated from the nascent 
neutron star. The ray-by-ray treatment essentially implies the 
construction of a ID transport solution in every lateral bin of 
the computational polar grid, assuming a spherically symmetric 
transport problem for the stellar conditions that are present in 
this angular bin. Therefore the ray-by-ray approach tends to un- 
derestimate the directional smearing of luminosity features that 
result from local and time-dependent emission increase in the 
semi-transparent accretion layer near the neutron star surface. 
Such features are, for example, caused by the SASI and convec- 
tion modulated mass inflow and must be expected to be more 
prominent in the plots showing the measurable luminosity L v> \, 
because this quantity is the isotropic equivalent luminosity re- 
sulting from neutrino emission occurring in one angular grid bin 
(see lefthand panels of Fig. |6]). In contrast, spatially localized 
effects will show up only with significant damping in the lumi- 
nosity L V 2, where neutrinos coming from the whole hemisphere 
facing the observer are added up with equal weighting (righthand 
panels of Fig. [6]). We consider these two evaluated cases as the 
extrema that define the theoretical limits of the true result. We 

4 The factor cos 6 accounts for the projection of the surface elements 
on the radiating sphere perpendicular to the line connecting source and 
observer. 
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suspect that L v \ yields the better approximation of the real situ- 
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Fig. 7. Mean energies of radiated neutrinos as functions of post-bounce time for our ID simulations (top) and 2D models (middle 
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that are averaged over the whole northern hemisphere of the computational grid. In all cases the evaluation has been performed in 
the laboratory frame at a distance of 400 km from the stellar center. 



magnitude of the luminosity fluctuations. It is important to note 
that the total energy loss of the supernova core in neutrinos (or 
one sort of neutrinos), which is usually defined as neutrino lu- 
minosity and computed by integrating the neutrino flux over all 
observer directions, shows SASI induced variability only on a 
much reduced level (see the corresponding plots in the Marek & 
Janka 2007 paper). This can be understood from the fact that the 
SASI fluctuations in the northern and southern hemispheres are 
phase shifted by roughly half a SASI period so that the emission 
maxima on one side fill the emission minima on the other. The 
combined signal is therefore much smoother than the neutrino 
flux that can be received by an observer in any chosen direction. 

The modulations of the apparent luminosity due to SASI and 
convective variations of the neutron star accretion rate in the 2D 
simulations are superimposed on the general luminosity evolu- 
tion that is characteristic of the employed progenitor star. These 
modulations add another feature on top of the effects that dis- 
criminate ID from 2D results and models based on the use of 
the L&S EoS from those with the H&W EoS. The 15 M s stel- 
lar model used in the present simulations shows a distinct drop 



in the mass accretion rate by roughly a factor of two when the 
steep density decline at the interface between the silicon shell 
and the oxygen-enriched Si layer falls through the stalled shock 
(see Appendices A and B of Buras et al. 2006b). This happens 
between 160 ms and 180 ms and leads to the pronounced de- 
crease of the luminosities of all neutrinos and antineutrinos at 
around this time. 

In Figs. [6] and Q we see an effect already mentioned in 
Sect. 13.11 namely that the proto-neutron star in the simulations 
with the softer L&S EoS is more compact and hotter and there- 
fore radiates higher luminosities and significantly higher mean 
energies of all kinds of neutrinos and antineutrinos (the mean 
neutrino energies (e v ) are defined here as the ratio of the energy 
flux to the number flux). This holds for ID as well as 2D models. 
Besides the quasi-periodic fluctuations, the 2D results exhibit the 
characteristic differences compared to the ID case that were dis- 
cussed before in much detail by Buras et al. (2006b): convective 
transport below the neutrinospheres leads to an accelerated loss 
of lepton number and energy, which manifests itself on the one 
hand in a slightly enhanced electron neutrino flux relative to the 
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electron antineutrino flux, and on the other hand in significantly 
increased muon and tau neutrino luminosities compared to the 
ID results. These differences are visible in both the polar lumi- 
nosities, L Vi ,i, as well as in the hemispheric ally integrated data, 
L Vi 2, of Fig. [6] On the other hand, the convective neutron star 
is less compact than its ID counterpart, for which reason the 
mean energies of the radiated neutrinos are slightly reduced in 
the 2D simulations (compare the upper two panels of Fig.|7]with 
the bottom ones). This phenomenon is particularly strong for the 
L&S EoS, in which case the more compact neutron star struc- 
ture seems to react more sensitively to the redistribution of lep- 
ton number and entropy associated with the convective activity 
in its interior. 

An interesting effect is the crossing of the mean energies of 
electron antineutrinos and heavy-lepton neutrinos ((ej; e ) becom- 
ing higher than (e v ,)) that occurs in five of the six panels of 
Fig. |7] at about 200 ms after bounce, but which cannot be ob- 
served until this post-bounce time for other energy moments; 
for example it is absent in the case of the rms energies of the 
energy spectra. These rms energies are ordered in the standard 
way at all times, with heavy-lepton neutrinos v x being more en- 
ergetic on average than v e , and the latter being more energetic 
than v e : (ey e ) rms < <ev e ) ms < (ey,> rms (this was pointed out 
already by Marek & Janka 2007). The energy crossing there- 
fore is a consequence of a subtle difference in the spectral shape 
of v e and v x . A close inspection of our results reveals that this 
effect is connected with the modification of the radiated neu- 
trino fluxes in the accretion layer around a successively more 
compact, contracting proto-neutron star. In contrast to the mean 
energies at a radius of 400 km, i.e. far outside of the neutrino 
source (see Fig. [7}, the corresponding data at a fixed density of 
10 11 gcrrT 3 always follow the relation (e^) % (e v ). This holds 
for our ID results and for the time-dependent unidirectional as 
well as (hemi)spherically averaged 2D results with both equa- 
tions of state; only during short phases of a strong increase of the 
accretion rate can the mean energy of electron antineutrinos ex- 
ceed the one of heavy-lepton neutrinos also at p — 10 11 gcm~ 3 . 
Similarly, accretion around the poles of the 2D models also 
causes the polar rms energies of v e and v x — measured out- 
side of the accretion layer — to converge ({ev e ) lms ^ ( 6 v I ) rms )> 
whereas the laterally averaged values are always well separated, 
in agreement with the hierarchy of the rms energies found in the 
corresponding ID simulations. Moreover, the neutrino radiating 
"naked" neutron star in spherically symmetric models with artif- 
ically initiated explosions does not produce the inversion of the 
mean v e and v x energies, whereas the accreting neutron star in a 
corresponding nonexploding model does (see Fig. 2 in Janka et 
al.2007). 

All these findings point to the importance of accretion for an 
explanation of the phenomenon. Accretion modifies the density 
profile outside of the neutrinosphere and thus is responsible for 
a subtle change in the v e spectrum relative to the v x spectrum. 
Two effects seem to contribute. On the one hand, with the inclu- 
sion of v x v A -pair production by nucleon-nucleonbremsstrahlung 
and by v e v e annihilation in state-of-the-art supernova models, the 
energysphere of heavy-lepton neutrinos is shifted to lower tem- 
peratures than in older simulations where these processes were 
ignored. The heavy-lepton neutrinos that leave their average 
neutrino-energysphere also experience a "spectral filter effect" 
in being down-graded in energy space by neutrino-nucleon and 
neutrino-electron collisions when diffusing through the cooler, 
outer layers of the neutron star (modern multi-energy trans- 
port calculations follow the corresponding spectral changes with 
great accuracy). For both reasons the spectra of v x and v e are 



now much more similar than in the older models (for a de- 
tailed discussion, see Raffelt 2001 and Keil et al. 2003). On the 
other hand, electron neutrinos and antineutrinos, which are emit- 
ted in rapid charged-current reactions of electrons and positrons 
with protons and neutrons, are more copiously produced in the 
hot accretion layer than muon and tau neutrinos. This leads to 
an enhancement of the v e flux at medium and high energies 
(^lOMeV) compared to the v x flux at these energies. In com- 
bination, these effects can raise the average energy of electron 
antineutrinos radiated from an accreting neutron star above that 
of heavy-lepton neutrinos. 

The described relative spectral changes become stronger at 
later times after core bounce so that a crossing of (e^) and (e Vx ) 
occurs. As the post-bounce time advances, the average ener- 
gies of muon and tau neutrinos grow slowly because the neu- 
trinospheric temperature increases in the shrinking neutron star. 
This increase is stronger for the soft L&S EoS than for the stiff 
H&W EoS (see Fig. |7). In contrast, the v e and v e created in 
the semitransparent accretion layer escape with mean energies 
that rise relatively faster, because accretion onto a more com- 
pact proto-neutron star releases more gravitational binding en- 
ergy. This leads to more extreme heating of the infalling and 
compressed gas and more abundant production of v e and v e with 
higher energies. This reasoning also explains why the polar-near 
data in Fig. [7] show the effect more clearly than the hemispher- 
ically averaged results: the SASI sloshing motions of the shock 
along the direction of the polar grid axis with alternating peri- 
ods of shock inflation and contraction lead to successive phases 
of gas expansion and mass accumulation on the one hand, and 
gas compression on the other, in particular around the north and 
south poles of the neutron star (compare Figs.[3]and[6]l. 

At first sight this argument seems in conflict with the obser- 
vation that the accretion funnels in the sequence of snapshots of 
Fig. |3] in particular the most prominent and persistent ones that 
carry most mass, reach down to the neutron star typically at inter- 
mediate latitudinal angles, i.e., around the equator; still the fluc- 
tuations of the neutrino emission are highest close to the poles. 
We understand this as a consequence of the following two facts. 
Firstly, because of the extremely large infall velocities only little 
neutrino loss occurs within the downdrafts. Instead, most of the 
binding energy is radiated when the accreted gas spreads around 
the neutron star and begins to settle onto the neutron star sur- 
face (for a discussion of this aspect, see also Scheck et al. 2006). 
Secondly, due to the SASI sloshing of the whole accretion layer, 
much of the gas accreted through the shock is redirected alter- 
natingly towards the north pole or the south pole. There the gas 
is decelerated near the grid (and sloshing) axis and deflected so 
that its duration of stay is longest in the polar regions. The gas 
that is sucked to a pole at an enhanced rate during the SASI half- 
cycle of shock expansion is compressed there during the next 
half-cycle of shock contraction. This is the reason why most of 
the losses of accretion energy happen around the poles. Future 
three-dimensional (3D) simulations with neutrino transport will 
have to show whether this behavior is a consequence of the as- 
sumption of axial symmetry or whether a similar effect can also 
be found in connection with the SASI in three dimensions. 

Concerning the possibility of measuring the SASI (and con- 
vection) induced luminosity variations by the detection of neu- 
trinos from a future Galactic supernova it is important to note 
that the strong compressional heating of accreted gas during the 
SASI half-cycles with shock contraction does not only boost the 
neutrino emission at the poles but also leads to the production 
of more energetic neutrinos. Therefore the variations of the neu- 
trino luminosities correlate and are in phase with modulations 
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Fig. 8. Fourier amplitudes of the time-dependent neutrino luminosities shown in Fig. [6] (the left/right column there corresponds to 
the left/right column in the present figure). In order to reduce the level of background noise, the evaluation was performed for post- 
bounce times t > 100 ms, for which the SASI modulation is clearly visible in the neutrino luminosities, and the signals were split 
in 17 overlapping segments of 128 ms, advanced in steps of 10 ms. The Fourier spectra of the ID models show a weak contribution 
from low-amplitude numerical fluctuations of the transport results. 



of the mean neutrino energies (see in particular the polar data in 
Figs.|6]and|7]l. While (for the polar emission) the luminosity dif- 
ferences between maxima and minima can be several 10% up to 
about 50% of the values in low-emission episodes and the mean 



energies show variations up to roughly 1 MeV or 10%, the ob- 
servationally relevant energy moments (approximately scaling 
like L (e) 2 or L (e), depending on the neutrino detector) fluctuate 
with an amplitude of typically ~50% and extrema up to nearly 
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100% of the values in the minima. The fluctuations are similarly 
strong for v e and v e , which are both abundantly created in the 
accretion layer, and somewhat reduced for heavy-lepton neutri- 
nos, whose emission is more indirectly affected by the pertur- 
bations of their decoupling region near the neutron star surface, 
e.g. through the surface gravity modes caused by the SASI ac- 
tivity and through changes in the density structure and optical 
depth associated with the mass motions in the SASI region. 

In Fig. [8] we provide the Fourier spectra of the neutrino lu- 
minosities displayed in Fig.|6l the lefthand panels belong to the 
polar luminosities, the righthand panels to the hemispherically 
averaged ones. The 2D simulations exhibit a broad peak be- 
tween about 20 Hz and roughly 400 Hz with a clear power ex- 
cess compared to ID results. Instead of the expected monotonic 
decline the Fourier transforms of the ID luminosities possess a 
flat shoulder at frequencies between ~50Hz and ~100Hz. This 
feature is a consequence of low-amplitude numerical noise pro- 
duced by the transport scheme, which is so small that it is physi- 
cally insignificant for the simulations. As expected from the time 
evolution of the shock radii and corresponding luminosity vari- 
ations, which both are larger in the case of the L&S run, the 
Fourier spectra of the 2D simulation with the L&S EoS dominate 
the result with the H&W EoS by a factor of 2-5 for all kinds of 
neutrinos. The flat local maxima in the ~25-80 Hz range, whose 
location roughly agrees with the peak of the Fourier amplitude 
of the t — 1 shock oscillation mode (Pig. |5j and also shows the 
tendency of slightly higher frequencies for the L&S case, are su- 
perimposed by a high-frequency component that originates from 
the faster convective modulations of the accretion flow and neu- 
trino emission. 



3.3. Gravitational waves 

Nonradial mass motions in the three regions of hydrodynamic 
instabilities as well as anisotropic neutrino emission are sources 
of gravitational-wave signals from stellar core-collapse events 
(e.g., Miiller et al. 2004 and references therein). We will analyze 
our data for EoS-specific properties of these signals now. 

The (quadrupole) gravitational-wave amplitudes, energy 
spectra, and spectrograms resulting from anisotropic mass mo- 
tion can be computed for axisymmetric models as described 
in Miiller & Janka (1997; Eqs. (10)-(12)), using the Einstein 
quadrupole formula in the numerically convenient form derived 
by Blanchet et al. (1990), and by standard FFT techniques. 
Assuming an observer that is located at an angle 9 with respect to 
the symmetry axis of the source, the dimensionless gravitational 
wave amplitude h{t) is related to the quadrupole wave amplitude 
Afo (measured in units of cm), the lowest-order non-vanishing 
term of a multipole expansion of the radiation field into pure- 
spin tensor harmonics (see Eq. (9) of Miiller 1997), according 
to 



time-dependent neutrino luminosity L v (t) of the supernova core 
and the emission anisotropy a v (f), 
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where R is the distance to the source. In the following we will 
always assume sin 2 9=1, i.e., the observer is positioned in the 
equatorial plane of the polar grid. The wave amplitudes due to 
the anisotropic neutrino emission are obtained with the formal- 
ism given in Miiller & Janka (1997; Eqs. (28)— (31)), which is 
based on the work of Epstein (1978). According to these ref- 
erences, the dimensionless gravitational-wave amplitude for the 
radiated energy in neutrinos can be written as convolution of the 
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(the quadrupole wave amplitude AS for neutrinos then results 
from inversion of Eq. ©), with 
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where dL v (Q, f)/dfi is the neutrino energy radiated at time t per 
unit of time and per unit of solid angle into direction dQ, and 
dQ. = -dcos&I^ defines the solid angle element in the polar 
coordinate system of the source. The angle-dependent function 
x ¥(9,(p) is given in Eq. (27) of Miiller & Janka (1997). For the 
case of an axially symmetric source, Kotake et al. (2007) derived 
a compact expression (see their Eq. (8) and note that our "F, fol- 
lowing the notation of Miiller & Janka 1997, is their cD), which 
is visualized as a function of the polar angle 9 in their Fig. 1 . For 
the further discussion it is important to note that has positive 
values in the polar cap regions between 9 = and 9 = 60° and 
between 9 = 120° and 9 = 180°, but adopts negative values in 
the equatorial belt for 60° < 9 < 120°. The combined polar caps 
and the equatorial belt have equally large surface areas. 

The quadrupole wave amplitudes of our two 2D models 
and the corresponding Fourier spectra are shown in Fig. [9] in the 
upper panels for the signals from aspherical matter flow, and in 
the lower panels for those from anisotropic neutrino emission. 

The matter signal consists of a superposition of quasi- 
periodic variations on timescales ranging from a few to several 
ten milliseconds. Two phases of enhanced activity can be distin- 
guished. The first one occurs within tens of milliseconds after 
core bounce and is followed by a more quiescent episode before 
a long-lasting period with an overall trend in growing amplitudes 
begins. The transient, early phase of strong gravitational- wave 
emission is a consequence of the prompt post-shock convection 
(cf. Sect. |3.1t and exhibits larger amplitudes, slightly higher fre- 
quencies (spectral peak at about 100 Hz instead of ~70Hz), and 
more power in our simulation with the stiffer H&W EoS (see 
Fig. [Toll. The reason for these differences is the fact that the re- 
gion of prompt convection involves more mass and extends to 
larger radii in the H&W case (Figs. [JJ and [2}, where also the 
growth rates of the instability are higher and the mass motions 
become more violent and faster. The prominent peaks of the 
Fourier spectra in Fig. [TOl correspond to a clearly dominant fre- 
quency that can be identified in the wave train during the first 
~80 ms after bounce. 

The later phase of intense gravitational-wave production is 
associated with the development of convective overturn and 
large-amplitude SASI oscillations in the post-shock layer as well 
as mass flows and convective activity in the proto-neutron star. 
Accordingly, the wave amplitudes at t 100 ms show variabil- 
ity with different frequencies. It is interesting to note that the 
enveloping curves of the wave amplitudes at this time follow 
roughly the evolution of the maximum shock radii in Fig. [4] 
This means that more active phases and stronger gravitational- 
wave emission correlate with larger shock radii. The nonra- 
dial gas motions in the post-shock region contribute mostly to 
the power radiated at frequencies up to about 200 Hz in accor- 
dance with the SASI Fourier spectra shown in Fig. [5] The main 
peak of the matter gravitational-wave spectrum at 300-600 Hz 
(H&W EoS) and 600-800 Hz (L&S EoS) as well as most of 
the lower "background" above / ~ 200 Hz, however, originates 
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Fig. 9. Gravitational-wave quadrupole amplitudes as functions of post-bounce time (left) and corresponding Fourier spectra 
{right) associated with mass motions (top) and anisotropic neutrino (v e + v e plus all heavy lepton neutrinos and antineutrinos) 
emission (bottom) for a distant observer in the equatorial plane of the axisymmetric source. Note that despite the trend toward a 
prolate shock deformation at late post-bounce times (see Fig. [3] and the lefthand panels of Fig. [5]l the neutrino amplitude develops 
negative values, suggesting that the signal is not determined by a dominant polar neutrino emission in this case. The prominent peaks 
of the matter signals at Fourier frequencies around 500 Hz and 700 Hz, respectively, are caused by mass motions in the proto-neutron 
star surface at densities between 10 11 g cirT 3 and 10 13 g cm 4 . 



from changes in the mass-quadrupole moment in deeper regions, 
namely at densities between 10 11 g cirT 3 and 10 13 g cirT 3 , where 
more mass is involved than in the post-shock layer. The main 
source of this signal are gravity modes and gas flows instigated 
by the violent impact of accretion funnels in the surface lay- 
ers of the nascent neutron star, which perturbs regions at den- 
sities up to p ~ 10 I2 gcirT 3 . A substantial part of the emission 
also comes from the convective zone inside the neutron star at 
p > 10 I2 gcm 3 . The importance of the interaction of the SASI 
and accretion flows with the proto-neutron star surface can be 
concluded from a closer inspection of the upper lefthand panel 
of Fig. [9]in comparison with Fig. Q] The gravitational-wave am- 
plitudes of the L&S run at t > 100 ms p.b. are significantly larger 
than those of the H&W model until near the end of our simula- 
tions. Only then the amplifying SASI and accretion activity be- 
gins to affect shells at p ~ 10 12 gcrrT 3 also in the H&W case, 
while earlier the involved regions in this model possess clearly 
lower densities. 



The maximum values of the dimensionless gravitational- 
wave strain h(f) for the matter signal are of the order of sev- 
eral 10~ 22 for a source at a distance of lOkpc, and the dominant 
peak may be marginally observable with LIGO I in the case of a 
Galactic supernova (for a detailed discussion of measurement as- 
pects, see Miiller et al. 2004). While the spectral peak is slightly 
higher if the nuclear EoS is soft and the proto-neutron star more 
compact, the stiffer H&W EoS leads to somewhat lower peak 
frequencies and thus to a shifting of the main emission closer 
to the region of highest sensitivity of the LIGO instrument. It is 
therefore not clear that a soft EoS is more favorable for a mea- 
surement of such gravitational-wave signals. 

Significantly larger amplitudes are generated by the 
anisotropic neutrino emission, however most power of this com- 
ponent of the gravitational-wave signal comes at frequencies be- 
low 100 Hz (Fig. [9}. The lower lefthand panel of Fig. [9] shows 
negative amplitudes with the continuous trend in growing ab- 
solute values and superimposed high-frequency modulations for 
both 2D simulations discussed in this paper. The negative val- 
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Fig. 10. Fourier spectra of the gravitational wave quadrupole am- 
plitudes of the mass motions caused by prompt post-shock con- 
vection during the first 100 ms after core bounce (see upper left- 
hand panel of Fig. |5). 



ues indicate a time-dependent quadrupole moment that is deter- 
mined by an excess of neutrino emission from regions near the 
equatorial plane, a fact that can be concluded easily from the lat- 
itudinal variation of the function ^¥(0) in the angular integral of 
the anisotropy parameter a v (f) of Eq. (0. This function becomes 
negative in the equatorial belt between 60° and 120°, see Fig. 1 
ofKotakeetal. (2007). 

A continuous decrease of the neutrino gravitational-wave 
amplitude to negative values is no common feature of all 2D sim- 
ulations. This behavior depends strongly on the particular emis- 
sion asymmetry that develops within the nascent neutron star and 
in its accretion layer. Miiller et al. (2004) found long-period vari- 
ations between positive and negative values, for example, for an 
explosion simulation of an 1 1.2 M Q star (cf. their Fig. 4) and for 
a 2D neutrino-cooling simulation of a convective proto-neutron 
star (cf. their Fig. 5), and they obtained a positive gravitational- 
wave amplitude with an overall trend in growth and only shorter 
periods of decrease for the neutrino emission of a centrifugally 
flattened, accreting neutron star at the center of a collapsing, ro- 
tating 15 M Q supernova progenitor (cf. their Fig. 3). 

This diversity of possible behaviors is in conflict with results 
obtained by Kotake et al. (2007), who reported almost monoton- 
ically rising, positive wave amplitudes for all of their 2D simula- 
tions of SASI unstable accretion shocks in collapsing supernova 
cores. They explained their finding by the particular properties 
of the SASI and the corresponding post-shock accretion flow in 
models that are constrained by the assumption of axial symme- 
try. In this case the SASI oscillations must occur along the polar 
grid axis, and Kotake et al. (2007) argued that the influence of 
this on the accretion flow makes the neutrino emission biased to- 
wards the regions around the symmetry axis (see also Sect. 13.21 ), 
thus leading to positive and growing wave amplitudes. 

Why do we come to different conclusions about the neutrino 
emission asymmetry than Kotake et al. (2007) in spite of our 
apparent agreement about the role of the SASI for the v e and 
v e production in the accretion layer? In the first place one must 
note that the modeling approach taken by Kotake et al. (2007) 
is fundamentally different from ours. Besides many other differ- 
ences, they, for example, did not simulate the collapse of a stellar 
progenitor model but started from a steady-state solution of an 



accretion shock with a mass accretion rate of 1 M Q s _1 ; they ex- 
cised the neutron star interior at a constant radius chosen as the 
inner boundary of the computational grid, where they assumed 
a fixed density of 10 n gcrrT 3 ; and they kept the fairly high 
mass accretion rate and the accretor mass constant during all 
of the simulation (time-independent boundary conditions at the 
outer grid boundary), using a spherically symmetric Newtonian 
gravitational potential and ignoring self-gravity of the matter 
in the accretion flow. Moreover and very important in the con- 
text of the present discussion of gravitational-wave signals as- 
sociated with anisotropic neutrino emission, they did not solve 
the energy-dependent, multi-flavor neutrino transport problem 
but employed a light-bulb description with prescribed, time- 
independent neutrino luminosities and spectra from the nascent 
neutron star and local source rates for lepton number and energy 
due to neutrino reactions in the accretion flow. 

Based on the possibilities available to us for comparison 
we cannot make a reasonable judgement whether any, and if 
so, which of these simplifications are causal for the discrepant 
gravitational-wave results obtained by us and by Kotake et al. 
(2007). However, we strongly suspect that the approach of the 
latter work is too simplistic for reliably estimating the emission 
asymmetry of neutrinos that is responsible for the gravitational- 
wave generation. This conclusion is supported by a close analy- 
sis of our neutrino transport data and the insight they yield into 
the formation of the neutrino emission asymmetry. Figure QT| 
provides the corresponding information for our 2D simulation 
with the L&S EoS. 

In the upper lefthand panel of Fig. QT| we display the time 
evolution of the neutrino gravitational-wave amplitude for this 
model. Different from Fig. [9] we have performed the evaluation 
now not only at a radius of 400 km, where the neutrino lumi- 
nosity has reached its asymptotic level, but also at the (time- 
dependent) radius where the density has a value of 10 11 gem -3 , 
i.e. in the vicinity of the electron neutrinosphere and thus below 
the accretion layer where neutrino cooling still enhances the out- 
going v e and v e fluxes. The upper righthand panel of Fig.fTTIcon- 
tains the corresponding parameters a v (t) of the neutrino emis- 
sion anisotropy, which exhibit rapid variations between positive 
and negative values on timescales of milliseconds to ^10 ms and 
growing amplitudes at later times. Both upper panels show an 
even clearer long-time trend to negative gravitational-wave am- 
plitudes for the neutrino luminosity leaving the neutrinospheric 
layer: a v (t) measured at 10 11 gcrrT 3 attains positive values only 
during the very short time intervals when the high-frequency 
modulations reach their largest amplitudes, while this param- 
eter at a radius of 400 km enters the positive side for sig- 
nificantly more numerous and longer time intervals, although 
C df'ff v (f') < for all times at both locations. Correspondingly, 

the gravitational-wave amplitude at p = 10 n gcnT 3 is even 
more negative than the one at r — 400 km, and their differ- 
ence (green line in the upper lefthand panel of Fig. fTTb sug- 
gests that the amplitude associated with the neutrino emission 
in the accretion layer is positive and reduces the neutrinospheric 
emission asymmetry. This seems to mean that the final emission 
anisotropy of the neutrinos is determined already by the flux that 
leaks out from the neutrinospheric region and in which the equa- 
torial region has obtained a slight flux excess. Below we will 
see, however, that this interpretation is not completely correct 
and the true situation is more complex. A comparion of the dif- 
ferent lines in the upper lefthand panel of Fig. QT| also reveals 
that the fluctuations of the radiated neutrino gravitational-wave 
signal with a period of 10-20 ms (frequency ~50-100Hz) orig- 




Fig. 11. Quantities that are relevant for the gravitational-wave signal associated with the anisotropic neutrino emission in our 2D 
model using the L&S EoS. Top left: Gravitational-wave amplitudes as functions of time after core bounce. The black curve was 
obtained by considering the neutrino radiation at a large distance from the source (identical with the corresponding line in the lower 
lefthand panel of Fig. [9] where, as in Figs. [6] and UJ the evaluation was performed in the laboratory frame at a radius of 400 km). The 
red line shows the wave amplitude computed from the neutrino luminosities that leave the neutrinospheric region (rather, a sphere 
where the density is 10 11 gcirT 3 ). The green line displays the difference between the black and red curves. Top right: Anisotropy 
parameter a v of the neutrino emission as function of post-bounce time. The color coding of the curves is the same as in the previous 
plot. Bottom left: Energy radiated in all neutrinos through the "polar cap regions" (E p ) and through the "equatorial belt" (E e ; see text 
for details), and the difference in these results (scaled by a factor of 10), integrated from the beginning of the 2D simulation shortly 
after bounce until the time f p b on the horizontal axis. As in the two panels before, the color coding of the curves corresponds to the 
two different locations where this analysis was performed. Bottom right: Total energy radiated in neutrinos through the "polar cap 
regions" (E p ), through the "equatorial belt" (E e ), difference in both scaled by a factor of 20, and sum of both (E* ot ) as functions of 
radius. The time integration was performed from the beginning of the 2D simulation shortly after bounce until the end of the model 
run at about 400 ms after bounce. 



inate from the accretion layer and are therefore connected to the 
S ASI and convective activity in the post-shock region and its in- 
fluence on the neutrino-cooling zone as described in Sect. 13.21 

3.3.1. Neutrino emission anisotropy 

In order to obtain deeper insight into the origin of the neutrino 
emission anisotropy, we have calculated the cumulative loss of 
energy in neutrinos (i.e., the time- and surface-integrated neu- 
trino flux) once for the polar cap regions and another time for 
the equatorial belt, again at r = 400 km and at p - 10 11 g crrT 3 . 
The results are shown in the lower lefthand panel of Fig. QT| 
The small differences E e - E p between equatorial and polar 



neutrino energy losses, scaled by a factor of 10, are plotted 
with dash-dotted lines. Again we see that at large distances as 
well as near the neutrinospheres the equatorial emission exceeds 
slightly the polar one during the whole evolution after t ~ 40 ms 
post bounce, and this excess is reduced by the neutrino pro- 
duction that occurs in the cooling layer outside of a density of 
10 n gcirr 3 . Separating the analysis for the different neutrino 
kinds, it becomes obvious that only the v e and v e luminosities 
grow in the cooling layer above the neutrinospheres, and that 
these neutrinos are more abundantly produced in the polar re- 
gions where the loss of accretion energy mainly occurs, as dis- 
cussed in Sect. 13.21 For electron neutrinos the emission in the 
accretion layer near the poles is so strong that their polar energy 
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loss at large distances can become bigger than the equatorial one. 
In contrast, the flux of heavy-lepton neutrinos, which are radi- 
ated from the neutrinospheric region with a persistent surplus 
of the equatorial energy loss, gets a bit reduced in the accretion 
and postshock layers between p = 10 11 g cirT 3 and r = 400 km, 
whereby their equatorial emission excess grows to roughly twice 
the value near the neutrinospheres. Their asymmetry determines 
the final emission asymmetry, because they contribute the ma- 
jor part (~56%) of the total neutrino energy loss (in fact, the 
pole-equator asymmetries of v e and v e at large distances have 
opposite signs and nearly cancel each other). These findings are 
in clear contradiction to the results obtained by Kotake et al. 
(2007), who reported that the neutrino gravitational-wave sig- 
nals in their models were dominated by the v e emission and the 
other neutrinos did only contribute very small amplitudes with 
the same signs (see their Fig. 3). 

What is the reason that slightly (2-3%) more energy is ra- 
diated from the neutrinospheric layer near the equator than near 
the poles? We find that the emergence of this asymmetry coin- 
cides with the development of the stable convective shell inside 
the neutron star. This might suggest that convective transport 
near the equatorial plane is slightly more effective than in the 
polar regions. We have no satisfactory explanation for this phe- 
nomenon, because one would expect convection to operate sim- 
ilarly in all directions and not to produce a global asymmetry 
on angular scales much larger than the size of the (nonstation- 
ary) convective cells. This, however, is a statement that is cor- 
rect in three dimensions but might not hold for convection in the 
2D environment of our simulations. Here convective plumes and 
downflows are toroidal features because of the axial symmetry, 
and the polar grid has a preferred direction along the polar axis. 
We suspect — without having performed a detailed analysis — 
that the operation of convective overturn with cells of such a 
toroidal structure near and along the polar axis might be slightly 
different from the in- and outflows near the equatorial plane with 
fluid motions perpendicular to the axis. If this were indeed the 
origin of the neutrino emission asymmetry, it would imply a se- 
rious warning that the computation of such anisotropies for es- 
timating gravitational-wave signals (and also pulsar kicks based 
on the recoil associated with anisotropic neutrino emission) re- 
quires 3D modeling, because in 2D the discussed effects may be 
just a consequence of the symmetry constraints. 

But is this fishy directional asymmetry really responsible 
for our prediction of the neutrino gravitational-wave signal? 
Fortunately, it is not! To see this, we have evaluated our results 
in yet another way. In the lower righthand panel of Fig. QT| we 
display the cumulative — i.e. time-integrated from the begin- 
ning until the end of our 2D simulation — neutrino energy loss 
(again summed up for neutrinos and antineutrinos of all flavors) 
as a function of radius for the polar, equatorial, and total emis- 
sion. The fourth (blue) curve represents the difference (E e - £ p ) 
between equatorial and polar terms, scaled by a factor of 20. 
Consistent with what we discussed above there is an excess 
of the equatorial losses until about 45 km. This radius encom- 
passes well the convection zone inside of the nascent neutron 
star during the whole post-bounce evolution of the 2D model 
with L&S EoS (Fig. Q]) and is larger than the location where 
p = 10 n gcirT 3 for post-bounce times t 170 ms (lefthand 
panel of Fig. @). At greater radii when one enters the accretion 
and SASI layers, however, the polar energy loss first dominates 
before the equatorial one wins again at r t 100 km. After a local 
maximum at a radius slightly larger than 100 km the equator-to- 
pole asymmetry settles to its final value. This radial variation is 
found to be qualitatively very similar for all kinds of neutrinos, 



although with different values of the local maxima and minima 
and of the final, asymptotic level at large radii. 

From the lower righthand panel of Fig.QT]we learn that the 
emission anisotropy carried by the outgoing neutrinos is not im- 
posed near the neutrinospheres and preserved outside, but it is 
completely determined by the accretion region behind the su- 
pernova shock and the SASI and convective inhomogeneities in 
this region. The final emission anisotropy just by chance gets 
close to the value near the neutrinospheres. Moreover, the sep- 
arate analysis for the different kinds of neutrinos reveals that it 
is the muon and tau neutrinos that govern the final anisotropy 
of the radiated energy flux. In the two models presented here 
the energy carried away by these neutrinos in the equatorial belt 
is slightly higher than the one escaping from the polar cap ar- 
eas. This as well as the total asymmetry, however, is likely to 
be a time- and model-dependent result and to depend strongly 
on the detailed structure of the inhomogeneities that evolve be- 
hind the deformed accretion shock. It could also depend on other 
model aspects such as the value of the mass accretion rate, the 
direction-dependent stagnation radius of the shock, or the rela- 
tive sizes of the fluxes of the different neutrinos. Such conjec- 
tures are not only supported by the H&W model in comparison 
to the L&S run: the neutrino gravitational-wave amplitude of 
the latter model exhibits a significantly clearer trend in a mono- 
tonic decrease to negative values, while in the H&W case one 
can see long phases where the amplitude or its time derivative 
become positive (Fig. |9). These conjectures are also in agree- 
ment with the results of other models described in the Miiller et 
al. (2004) paper, and, last but not least, they are suggested by 
the fact that the quantity E e - E p in the lower righthand panel of 
Fig- Ellis a quantity with large radial variations in both positive 
and negative directions. It is rather unlikely that such a complex 
radial behavior can be generic for models in two dimensions so 
that in all cases and at all times the same effects are recovered 
that ultimately determine the emission asymmetry in our mod- 
els, namely (i) a delicate almost-cancellation of the combined 
asymmetry of v e and v e due to the production and absorption 
processes of these neutrinos in the neutrino-cooling and heat- 
ing regions, and (ii) a final dominance of the anisotropy of the 
radiated muon and tau neutrinos as a consequence of their asym- 
metric transport through the semi-transparent accretion layer, in 
which the total heavy-lepton neutrino luminosity hardly changes 
but the equator-to-pole asymmetry grows. 

Our analysis therefore reveals that the gravitational-wave re- 
sults for neutrinos depend on physical effects that are either 
not included or treated in a much simplified manner in previ- 
ous works. It further reveals that a reliable and quantitatively 
meaningful computation of the tiny neutrino emission asymme- 
try is extremely difficult and does not only require a full transport 
treatment and good numerical resolution in all relevant regions, 
but ultimately must be done on the basis of 3D models. 

4. Summary and conclusions 

We have presented results of two stellar core-collapse simula- 
tions performed with the Prometheus- Vertex code for a non- 
rotating 15 M Q progenitor star in axial symmetry. One of these 
models was computed with a soft version of the Lattimer & 
Swesty (1991) nuclear equation of state, and the other one with 
the stiffer Hillebrandt & Wolff (1985) EoS, which produces a 
significantly less compact proto-neutron star. For these simula- 
tions we have compared the neutrino and gravitational-wave sig- 
nals between core bounce and 400 ms later. In particular we were 
interested in signal features that may yield information about the 



18 



Marek et al.: Equation-of-state dependent supernova signals 



action of the standing accretion shock instability (SAS1) in the 
supernova core during the long phase of post-bounce accretion 
that precedes the onset of the supernova explosion in all recent 
multi-dimensional calculations of stellar core collapse. 
Our main results can be summarized as follows: 

(1) In both simulations we found that the weakened bounce 
shock produces a negative entropy gradient that partly over- 
laps with the negative lepton number gradient left behind by 
the escaping shock breakout burst of electron neutrinos. This 
region is Ledoux unstable and develops prompt post-shock 
convection, which leads to an early burst of gravitational- 
wave emission lasting several ten milliseconds after core 
bounce. Its frequency spectrum shows a prominent peak 
around 70-100 Hz. 

(2) A later, long-lasting phase of strong gravitational-wave 
emission with a long-time trend in growing wave amplitudes 
sets in at t ^ 100 ms after bounce. It is produced by nonra- 
dial mass motions and anisotropic neutrino emission that are 
caused by convective overturn inside the proto-neutron star 
and in the neutrino-heating layer behind the stalled shock, 
and in particular by vigorous sloshing of the accretion shock 
due to the development of low-^ mode SASI activity. The 
wave train for the matter signal consists of a superposition of 
quasi-periodic variations on timescales from a few millisec- 
onds to several ten milliseconds. The spectra of the dom- 
inant I = 1,2 SASI modes peak between about 20 and 
roughly 100 Hz, and the gravitational-wave emission also 
shows significant power at frequencies up to ~200 Hz. The 
most prominent maxima of the gravitational-wave spectra, 
however, are located at higher frequencies of 300-800 Hz 
and originate mainly from the outer layers of the neutron star 
where rapid gas flows are triggered by the supersonic impact 
and deceleration of accretion funnels. 

(3) The high-frequency peak of the gravitational-wave spec- 
trum and less clearly also the secondary maximum near 
the SASI frequency exhibit a dependence on the compact- 
ness of the nascent neutron star and thus on the proper- 
ties of the high-density EoS: the more compact neutron 
star leads to more powerful shock oscillations and signifi- 
cantly larger gravitational-wave amplitudes earlier after core 
bounce. Also the characteristic frequencies of the low-^ 
SASI modes and of the gravitational-wave signal are higher. 
While the main peak of the wave spectrum is located at 300- 
600 Hz for the stiff er EoS, it can be found at 600-800 Hz for 
the softer EoS with the more compact remnant. 

(4) The more compact neutron star also radiates neutrinos with 
higher luminosities and greater mean energies. Anisotropic 
neutrino emission produces a low-frequency component of 
the gravitational- wave signal at / S= 100-200 Hz with an 
amplitude that dominates the matter signal by up to a fac- 
tor of two. The emission anisotropy is mostly established in 
the cooling and heating regions of the accretion layer outside 
of the neutrinospheres and is determined by the anisotropic 
transport of muon and tau neutrinos through this layer, which 
leads to a slightly higher flux of muon and tau neutrinos near 
the equatorial plane of the polar grid, i.e., perpendicular to 
the direction of the SASI oscillations of the supernova shock. 
Despite the stronger emission of v e and v e from the accre- 
tion layer in the polar regions (i.e. in the direction of the 
SASI shock expansion and contraction), the gravitational- 
wave amplitude associated with the anisotropic energy loss 
in neutrinos exhibits a nearly monotonic trend to negative 
values. 



(5) The neutrino luminosities and mean energies, in particular 
the ones that can be received by an observer in the polar (and 
thus SASI) direction, show a quasi-periodic time variability 
with an amplitude of several ten percent (up to about 50%) 
of the minimum values for the luminosities and of roughly 
1 MeV (up to 10%) for the mean energies. The luminosity 
fluctuations are correlated and in phase with the energy vari- 
ations and are somewhat bigger for v e and v e than for heavy- 
lepton neutrinos. They originate from alternating phases of 
gas accumulation and compression in the accretion layer in 
particular around the poles, caused by the expansion and 
contraction of the shock in the course of the SASI oscil- 
lations. The frequency spectra of the neutrino luminosities 
show most power between about 20 and 200 Hz and signif- 
icant power in a decaying tail up to about 400 Hz as a con- 
sequence of emission variability associated with faster, con- 
vective modulations of the accretion flow between supernova 
shock and nascent neutron star. Since SASI and convective 
activity occur inseparably in the gain layer and stir the cool- 
ing layer below, power is found to be distributed over a wide 
range of frequencies in the Fourier spectra of the luminosi- 
ties. 

The SASI induced temporal modulations of the luminosities 
and mean energies are features by which neutrinos can provide 
direct evidence of the dynamical processes that occur in the su- 
pernova core. Only few other cases are known where neutrinos 
carry such information, for example the prompt flash of electron 
neutrinos as a signature of the shock breakout from the neutri- 
nosphere, or a possible second neutrino burst or abrupt termina- 
tion of the neutrino emission, which are expected if the neutron 
star collapses to a more compact object or a black hole. A de- 
tection of the predicted characteristic, quasi-periodic variation 
in the neutrino event rate from a Galactic supernova would con- 
firm the action of the SASI in the supernova core and would thus 
provide extremely valuable support of our present understanding 
of the core-collapse physics: the SASI is discussed as possible 
origin of global supernova asymmetries and pulsar kicks and as 
a potentially crucial agent on the route to the explosion. 

Moreover, nonradial mass motions such as convection 
and the SASI in the supernova core after bounce produce 
gravitational-wave emission with significant intensity. We found 
that the wave spectrum exhibits a broad hump between about 
20 and 200 Hz and a main maximum with clearly more power 
at considerably higher frequencies. The latter is associated with 
SASI and accretion induced gas flows in the neutron star surface 
layers, where much more mass is involved than in the fairly di- 
lute post-shock region. Both the primary and the secondary max- 
imum of the spectrum are sensitive to the structure of the neutron 
star: larger gravitational-wave amplitudes and correspondingly 
higher spectral maxima as well as higher frequencies are ob- 
tained with a softer equation of state and more compact remnant. 
The relative SASI variations of the neutrino luminosities do not 
appear to be very sensitive to the nuclear equation of state, and 
a clear EoS-dependent difference in the modulation frequencies 
is also not visible, but the absolute values of the luminosities 
and mean energies of the neutrinos radiated by the more com- 
pact and hotter, accreting proto-neutron star in our simulations 
are 10-20% higher. 

Anisotropic neutrino emission is also an important source of 
gravitational waves with amplitudes that can dominate the mat- 
ter signal. The corresponding spectrum possesses most power 
at frequencies below 100-200 Hz. The emission anisotropy is a 
function with rapid short-time variability (with periods of mil- 
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liseconds to ^;10ms) and with an overlying, nearly monotonic 
long-time trend toward slightly stronger emission from equato- 
rial regions, i.e., perpendicular to the direction of the S ASI shock 
oscillations. This trend decides about the sign of the neutrino 
gravitational-wave amplitude. A detailed analysis of our mod- 
els reveals that it is the result of complex transport physics in 
which the neutrinospheric flux, which has a low pole-to-equator 
anisotropy, is modified strongly in the inhomogeneous environ- 
ment of the accretion layer around the forming neutron star. The 
pole-to-equator asymmetry changes with radius several times in 
both directions and the contributions from v e , v e , and heavy- 
lepton neutrinos partly cancel each other. A reasonable determi- 
nation of these transport asymmetries, which is not only impor- 
tant for gravitational-wave calculations but also for pulsar kick 
estimates based on neutrino recoil, requires detailed transport 
calculations. Even our sophisticated ray-by-ray approach may 
not be sufficient for accurate predictions of measurable effects, 
but a full multi-angle treatment of the transport may be neces- 
sary. 

From this discussion it is clear that any reliable assess- 
ment of the detectability of SASI effects in the neutrino and 
gravitational-wave emission will ultimately require results from 
three-dimensional models. The growth of SASI modes with non- 
vanishing azimuthal order m was shown to be of potential impor- 
tance, and the presence of such nonaxisymmetric asymmetries 
can have a big influence on the structure and the dynamics of 
the accretion layer (Blondin & Mezzacappa 2007, Yamasaki & 
Foglizzo 2008, Iwakami et al. 2008). It is probable that also the 
emission asymmetries and time modulations of the emitted sig- 
nals exhibit significant differences compared to 2D models, in 
particular because in 3D the SASI axis is not forced to coincide 
with an axis of symmetry and with the axis of the numerical grid. 

The SASI and convective activity, its duration as well as 
strength, must be expected to vary with the progenitor star and 
to depend on the length of the postbounce accretion phase un- 
til the onset of the explosion (for hints to that, see the results 
for 11.2M and 15 M Q explosions in Marek & Janka 2007). 
In the nonlinear stage, the SASI and convective motions are of 
chaotic nature and therefore stochastic differences also in the 
large-scale pattern and global asymmetry are likely to develop 
from small initial differences even in progenitor stars (models) 
of the same kind (see, e.g., the different explosion asymmetries 
obtained in a large set of 2D simulations performed by Scheck 
et al. 2006). Nevertheless, basic features of the instabilites, in 
particular their mere presence and long-time growth behavior 
and the sizable modulation they impose on the emitted signals 
over a period of possibly hundreds of milliseconds, should be 
common to different realizations. Judged from our present un- 
derstanding, however, a template-like uniformity of the signal 
structure with a well defined frequency (or frequency range) and 
amplitude that exhibit a regular dependence on the explosion and 
progenitor properties, appears not to be very likely. More model 
calculations are therefore needed to find out which signal prop- 
erties such as the modulation frequency, are sufficiently generic 
to yield information about the neutron star equation of state. 

4.1. Detectability of SASI modulations 

Should the SASI signatures and their magnitude discussed in this 
work be confirmed by 3D results, the question arises whether 
they will be measurable with existing or emerging neutrino and 
gravitational-wave experiments in the case of the next Galactic 
supernova. A detailed prediction of supernova neutrino signals 
in terrestrial detectors is a highly complex problem, because 



the neutrinos radiated by the nascent neutron star undergo fla- 
vor conversions so that the flavor arriving at the Earth is not 
necessarily the same as the one leaving the source. On the one 
hand, the flavor evolution is affected by neutrino-matter oscilla- 
tions during the neutrino propagation through the density gradi- 
ents in the supernova mantle and envelope (where the H and L 
resonances are encountered at densities around 10 3 gcrrT 3 and 
10 gem 3 , respectively) and through the Earth. On the other 
hand, neutrino conversion is also induced by the nonlinear self- 
interaction associated with neutrino-neutrino forward scattering 
in the dense neutrino background close to the neutron star (for 
an application of these phenomena in the context of the neu- 
trino signals from oxygen-neon-magnesium core supernovae, 
see Lunardini et al. 2008, and from black hole formation in stel- 
lar collapse events, see Nakazato et al. 2008; for a discussion 
of the flavor conversion physics, see the references to original 
publications in the latter two papers). The flavor composition of 
the neutrino burst in the detector and thus the detector response 
therefore do not only depend on the type and structure of the 
progenitor star but also on a number of degrees of freedom as- 
sociated with unknown neutrino properties and the location of 
the experiment, e.g. the neutrino-mass hierarchy (normal or in- 
verted), the still undetermined mixing angle between the mass 
eigenstates 1 and 3, and the nadir angle of the neutrino path 
through the Earth. 

Because of the involved complexity we defer an investiga- 
tion of the flavor conversion effects and their consequences for 
the neutrino measurements to a future paper. However, since the 
SASI modulations affect the luminosities and mean energies of 
neutrinos and antineutrinos of all flavors with similar strength, 
we can ignore the flavor oscillation effects for a simple esti- 
mate to convince us that the SASI modulations of the radiated 
signal should be detectable in principle. Super-Kamiokande, for 
example, is expected to capture roughly 8000 electron antineu- 
trinos through their absorption by protons, if a Galactic super- 
nova occurs at a distance of lOkpc and releases at total energy 
of 5 x 10 52 erg in v e with an average energy of 15 MeV. A lumi- 
nosity of 4 x 10 52 erg s" 1 , which is a typical value for the emis- 
sion maxima during the SASI phase, should account for an event 
rate of approximately 7000 v e captures per second, correspond- 
ing to about 70 events in the ~10ms intervals of peak emission. 
If L v (e 2 )(e v ) _1 is about 50% lower in luminosity minima, about 
one half of this event rate and event number can be expected 
during the low-emission periods. Therefore the SASI modula- 
tions will lead to a variability of the measured signal that can be 
larger, though not much, than the statistical -^V v> d e t fluctuations, 
iVy.det being the number of detected neutrinos in a time interval 
of 10 ms. For Hyper-Kamiokande, whose fiducial volume would 
be roughly 30 times bigger than that of Super-Kamiokande, the 
Poisson noise would be reduced to a few percent of the event 
number within 10 ms, so that the SASI modulations should cre- 
ate a clear signature. Even more powerful for this purpose is 
the IceCube experiment. It will be able to identify the supernova 
neutrino burst by the Cherenkov glow caused by the neutrino en- 
ergy deposition in the ice, which will show up as time-correlated 
noise in all phototubes. For a supernova at a distance of lOkpc, 
the v e luminosity maxima during the SASI phase are expected 
to cause a photon count rate of more than 10 6 per second in the 
4800 optical modules (see Dighe et al. 2003). The -1.2 x 10 4 
events within 10 ms intervals and the variations of several 1000 
events between periods of luminosity peaks and minima (in the 
case of IceCube these variations scale with L v (e 3 )(e v ) _1 , because 
the crucial quantity is the energy deposition in the detector) dom- 
inate by far the Poisson fluctuations associated with the back- 
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ground counting rate of 300 Hz per optical module; the Poisson 
variability of the background corresponds to JW^ ~ 120 photon 
events in all phototubes, which is ~1% of the supernova signal in 
a 10 ms interval (an impression of the detection capability even 
for a supernova model with a lower fluctuation amplitude of the 
neutrino emission than in the models discussed here, can be ob- 
tained from the upper lefthand panel of Fig. 2 in Dighe et al. 
2003). 

The detectability of gravitational-wave signals very similar 
to the ones presented in this work was investigated in much detail 
by Miiller et al. (2004). Due to the restriction to a 90-degree grid 
(one hemisphere) in the case of the 15 M core-collapse model 
sl5r considered in the latter paper, the contributions of dipolar 
SASI modes were not included in there. Moreover, since the col- 
lapse evolution was followed to later times in our most recent 
simulations (about 400 ms instead of ~270 ms after bounce) and 
the quadrupole amplitudes exhibit a trend in rise during the later 
post-bounce phases, the signal-to-noise (S/N) ratios and maxi- 
mum detection distances given for LIGO I and Advanced LIGO 
by Miiller et al. (2004) are likely to be only lower limits. Since 
the S/N ratio scales roughly linearly with the Fourier transform 
of the wave amplitude, the detection prospects grow with longer 
integration periods of the signal and thus higher Fourier power. 
The measurement of dominant gravitational-wave power at fre- 
quencies of several hundred Hz with growing strength over a pe- 
riod of hundreds of milliseconds would therefore be a strong in- 
dication of SASI and convective activity around the nascent neu- 
tron star. Miiller et al. (2004) estimated a S/N ratio of roughly 70 
in Advanced LIGO for a 15 M Q supernova at 10 kpc, correspond- 
ing to the possibility of capturing gravitational waves from such 
sources up to a distance of roughly 100 kpc. A third-generation 
underground interferometer facility (Einstein Telescope) would 
lead to a gain in sensitivity of approximately a factor of 20 in 
the most relevant frequency range between 100 Hz and 1000 Hz 
(Spiering 2007). This would increase the expected detection rate 
of supernovae by up to a factor of 3 (Ando et al. 2005). A de- 
tailed analysis of the gravitational-wave signals and their ob- 
servational implications from a greater sample of core-collapse 
models will be presented in a forthcoming paper. 
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